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BE INTRODUCTION 

The enhancement of synoptic-scale precipitation by topography is well understood, 
and has been routinely verified. Fine-mesh mesoscale models have proven particularly 
accurate in reflecting orographic influences on large-scale precipitation events (Colton 
1976, Collier 1975). Simpler estimates of this enhancement factor have been shown to 
be almost exclusively functions of regional topography combined with low-level wind 
fields. synoptic-scale vertical motion. and a knowledge of the low-level relative humidity 
(Browning et al. 1975, HUl et al. 1981). 

Thus far, only detailed mesoscale models have been able to reproduce the small- 
scale orographic enhancement process. With synoptic-scale precipitation products 
showing little or no regional structure, and mesoscale forcing fields generally unavailable, 
forecasters rely mostly on local knowledge and vague statistical multipliers to predict 
areas of heavy rainfall. Consequently, the ability to quantify local orographic effects and 
their associated distribution of precipitation based solely on synoptically derived wind and 
moisture fields would be of particular use to the forecaster. 

In this study, the above mentioned quantities of subsynoptic-scale topography and 
synoptic-scale wind and moisture data are combined in an effort to reproduce areas of 
enhanced precipitation during a series of eastern Pacific extratropical cyclones which 
passed over California during the period 1200 UTC 27 February 1991 to 1200 UTC 04 


March 1991. Topography. conmibimesl sth ereh lo erre pre fre 1.1 us CIEN by the 


National Meteorological Center's (NMC) final analyses, produces a distribution of the 
total vertical moisture flux. This vertical moisture flux estimate is correlated with hourly 
precipitation data collected by the National Oceanic and Atmospheric Administration 
(NOAA), covering 244 reporting stations throughout California for the same period, to 
produce an equivalent rainfall amount. 

The overall intent of this study is to aid forecasters in identifying areas most likely 
to experience orographic enbancement of precipitation for a given synoptic event and 
provide a quantitative precipitation estunate. This procedure may be of use in estimating 
seasonal snowpack levels in otherwise inaccessible regions. as well as providing near term 
flood potential for specific drainage basins, while only requiring synoptic-scale 
information. More specific objectives of this study are to: 

- Assess the degree of accuracy to which precipitation distribution can be represented 
using mesoscale topography and large-scale wind and moisture fields. 


- Develop a simple quantity, calculated locally from model output products, to 
represent the enhancement of precipitation due to topography. 


- Present specific synoptic-scale conditions and/or regional topographic characteristics 
which either promote or discourage the use of this parameter. 


- Determine the level of improvement, if any, the employment of this parameter 
shows over current precipitation forecast products. 
Following a discussion of background work done in the field of orographic 
precipitation enhancement, an outline of the dynamical basis for the vertical moisture flux 
as a precipitation parameter is presented, az we Hors tetef de coy cer af the required 


input fields. Next, this parameter is correlated with observed precipitation ta ebtain a 


statistical fit assigning specific precipitation accumulations to the calculated values 
followed by an examination of the factors which affect correlation. Results obtained by 
applying the precipitation parameter both generally and regionally are examined, including 
a comparison to current precipitation forecast aids (specifically, the Nested Gnd Model 
(NGM) 12 hour precipitation forecast). Conclusions focus on the usefulness of this 


technique as a forecast aid. 


J. BACKGROUND 

It is certainly no secret that mountainous areas often receive greater amounts of 
precipitation than do upwind or downwind low-lying regions. However, it was not until 
Bergeron (1965) introduced his "seeder cloud - feeder cloud" concept that a specific 
mechanism for topographically enhanced precipitation became accepted. As described by 
Raddatz and Khandekar (1979) this process consists of low-level (feeder) clouds forming 
in upslope areas which "wash out" when seeded by raindrops falling from synoptic-scale 
(seeder) clouds above. The idea is that the feeder cloud is constantly replenished and will 
continue to enhance precipitation in proportion to its depth, provided seeding from above 
js uninterrupted. Storebo (1976) and Bader and Roach (1977) were able to confirm 
Bergeron’s theories by results which related a 10% increase in precipitation to a 1000 m 
feeder cloud depth. Interestingly, Browning et a/. (1974) suggested that in some instances 
(namely, cases of mid-level potential instability), orography alone could be responsible 
for generating precipitation, even if synoptic-scale seeding was not available. According 
to Browning et al. (1975), condensation over mountainous regions is roughly proportional 
to the low-level wind velocity component normal to the topographic barrier combined 
with synoptic-scale vertical motion and moisture quantities. This relationship was 
generally verified by the mesoscale numerical model developed by Collier (1975). 

As mentioned in the introduction, Colton (1976) achieved remarkably accurate 


results using a fine-mesh (4.3 KM grid) model, o voto lar o LO 


to tbe Sierra Nevada of California. A key assumption made by Colton was to define the 
orographic precipitation as a quantity dependent on the presence of synoptic-scale 
precipitation. As such, his scheme is only able to enbance synoptic-scale precipitation, 
and cannot force precipitation based solely on orography. 

A. later study by Hill et al. (1981) further confirmed the seeder-feeder enhancement 
theory. They additionally showed that if the total precipitation rate exceeded a given 
threshold value (0.5 mm hr’), the topographic enhancement factor was not dependent on 
tbe background precipitation rate. In short, as long as the background rainfall rate was 
sufficiently high, enhancement was simply proportional to the low-level moist flow 
interacting with topography, and was not a function of the amount of seeding from above. 
This carries the important implication that orographically induced precipitation can simply 
be added to synoptic-scale rainfall to obtain a total, rather than the orographic term being 
a function of the synoptic-scale tenm. This relationship provides the basis for deriving 
the precipitation parameter described in the following chapter. 

More recent studies have attempted to modify feeder-seeder theory to account for 
additional effects including stratification and wind drift. While showing the effect of 
Stratification to be negligible as it relates to total washout, Carruthers and Choularton 
(1983) did show wind drift effects, which were not accounted for in models such as Bader 
and Roach (1977), to be detrimental to both local maxima and total enhancement by 
"advecting" orographic precipitation to downwind locations. Additionally, this study 
MAS AT DE 1deagpieseniestiotse s Ense E a a e a ndependent of 


seeding rate (also contradictorv to the Bader and Poach (1977y model) 


Most recently, Stow, Bradley and Gray (1991) conducted a study with results 
largely consistent with seeder-feeder theories, and emphasized the significant differences 
between the previous enhancement studies. They regarded the work so far conducted as 
inconsistent in terms of dealing with the effects of the low-level wind speed, seeding rate, 
and barrier height. Moreover, they felt that although the basic physics of topographic 
enhancement are well understood, refinement is required to handle the above 
inconsistencies. The present study js not concemed with any of these "refinements", but 
rather with using the existing theory to derive a simplified method of representing 
mesoscale-resolution topographically enhanced precipitation from exclusively synoptic- 
scale data. 

It should be noted that topography plays a larger role in weather than simply 
enhancement of precipitation. In fact it has been shown that topography can influence 
the low- to mid-level flow in such a manner as to completely reposition large-scale 
pressure centers (Sawyer 1959), as well as greatly altering the small-scale wind patterns 
(Danard 1977). An attempt to combine small-scale deformation of the wind field due to 
topography with established enhancement theory was undertaken by Anderson and Nilsson 
(1990) and had only moderate success, mostly owing to large errors in obtaining 
precipitation measurements in drifting snow. Any deformation of the wind field by 
topography will be ignored in this thesis, in terms of calculating our precipitation 
parameter. As such, a main objective of this study is to determine the feasibility of 
representing mesoscale precipitation Un CI OO TG Mc ac N 


mesoscale wind fields. 


JU PROCEDURE 


A. “THEORY 

Preceding studies have estimated orographic enhancement of precipitation by 
determining the vertical extent of a feeder cloud which forms in upslope regions based 
on low-level winds and wet bulb potential temperature. This cloud depth can be directly 
reJated to orographicallv induced precipitation provided seeding from above (synoptic- 
scale rain) is available. Assuming the “washout” efficiency is properly handled, numerical 
models can easily determine the enhancement factor based on calculated feeder cloud 
depth (Raddatz and Khandekar 1979, Bader and Roach 1977). 

Since the purpose of this study is to provide an easily-calculable means of 
estumating topographical enhancement of precipitation, a slightly different approach than 
a numerical model simulation is taken. This estimate of topographic precipitation 
enhancement is based on a "ground-up" quantity reflecting the total saturated moisture 
flux at a given height (assumed to be the cloud base), (wq),. Conceptual restraints are 
imposed such that any upward saturated moisture flux at this level is condensed, and not 
allowed to either be redistributed into higher levels of the atmosphere, or result in 
supersaturation. Essentially, any such flux is completely "washed-out" as precipitation. 
Since this flux is representative of precipitation only under saturated conditions, a 
e 700 mb relative bumidity volue ef 755. € S ocploetee needy ts used to 


represent saturation at the svnoptic scale) will be aprbied te the fw, field. 


B. DERIVATION OF PRECIPITATION PARAMETER 
Bader and Roach (1977) used the following form of the continuity equation to 


represent the liquid water content: 


P41 -C-A 
Dt (1) 





where Dq,/Dt represents rate of change of liquid water content following the air parcel 


trajectory. C is used to represent the condensation rate while A is the rate of washout 


from above. Jt follows that Dq /Dt 1s simply the total precipitation minus evaporation 


under saturated conditions, such that 


e (EN) 
Da 


(2) 


where P and E stand for precipitation and evaporation, respectively. Following Holton 


(1979, pp 40-44), 
dq, 0q, | dg. 9. dq, 
= u +y +W 


C 


de or | ox dy op 
represents the total time derivative of g, in terms of the local tendency and advective 


components. Rearrangement in vector notation and use of the isobaric form of the 


continuity equation with (2), gives 


oq, 
a +V-(g,0)=-(P-E). 
dt 


Changes in the saturation specific humidity (q,) result from temperature and pressure 


changes. which are largest due to vertical advection. Thus, the local tendency dq /ot is 


smaller than the flux divergence V-gU and will be ignored in this calculation. The P-E 
term is the net condensation, which 1s equivalent to the rainfall when integrated over the 
depth of the atmosphere. The net condensation is written simply P in the following 


equations. We are left with 


| (4) 
=P=V 190), 


as a first-order approximation to the net precipitation. Separating (4) into horizontal and 


vertical components yields 








qu) , 0(qv) (qo)‏ م_ 
ox dy op‏ 


This can be integrated from the surface to the top of the atmosphere H, to give 





te / dqu) dqy) (5) 
- [Pap +(qw), =(90), 4 EB tp 
Po ١ j Po ox dy 


Py 
where (qw), represents total moisture flux at H, which is zero. The term - [Pap 


Po 


represents the rainfall. R, and the right-hand side of equation (5) gives the contribution 
Pr 


from both topographic uplift, (q0), , and synopUc-scale forcing, + [9 QU)]ap. The 


Po 


synoptic-scale forcing is simply the net horizontal flux divergence of saturated air into the 


column. Because q, is largest at low levels, the net moisture flux is dominated by the 


lower part of the atmosphere. Consequently, this term is approximated by the integrated 
moisture flux convergence from the surface to 700 mb. The approximate rainfall is then 


given by 
700 
R AR. flv Xd Op, 
Po 


with R representing net precipitation in the column Roth terms appearing on the right- 


hand side of the equation are easily obtained hero standard wind am] moisture data 


10 


combined with surface topography. The 40 km resolution topography used in this study 


dh dh 


is shown in Figure 1. Specifically, the term (90), is given as q, Ces w, 
0 0 y 


0x 


with u, v and h representing horizontal wind velocity components and surface slope, 


respectively, while q, refers to the surface moisture value. 
0 


The surface forcing (qw). , and ulumately (qq), , was calculated using the 1000 mb 


wind field and topography in meters. with the 1000 mb winds treated as crudely 


representative of synoptic-scale surface winds for this purpose. The resultant surface- 


forced moisture flux is in units of ¢/kg-cm/s. An approximate conversion to pressure 


coordinates is made by applying the relation qw=-pgwyg. 


TOU 
The synoptic scale term )+ | 1990014 is calculated using a centered finite 


Fo 


difference scheme to compute the moisture flux divergence at levels specified by the wind 
and moisture field inputs. Summations over the desired vertical distance yield the net 


moisture flux divergence in the surface to 700 mb layer, which is the vertical moisture 


flux g@ through the 700 mb level. This is the quantity w0f, referred to in the Results 


chapter. Figures 2 and 3 show sample contoured surface vertical velocity and qw flux 
ME A Obtained HOM ece mme hoda. Uli topes iii repase zenn an Ergine 2 has 


been subjected to a five passes through a three-peumt smoothing process to ensure thus 


11 


quantity is of the same order of magnitude as the synoptically derived forcing in Figure 


3. Finally, the total moisture flux at 700 mb, (go),,, (hereafter referred to as wqf), is 
obtained by simply adding the topographic and synoptic-scale components, utilizing the 


rough equivalency of p bars! and cms’! as described above (Figure 4). Thus both w0f and 


wgf are in units of ekg "1 ubars "1. This unit definition applies to all figures showing wOf 
or wqf fields unless otherwise noted. 


Note that this (q9),,, quantity is not in units of rainfall. The actual associated 


rainfall is obtained by frst employing a 75% 700 mb relative humidity threshold to the (@q),,,, 


field to ensure saturation. This field is then correlated with observed 12 hour 


accumulated precipitation centered on the analysis time. A linear regression is calculated, 


yielding a conversion factor based on the slope and y-intercept, allowing(q),,,, 


quantities to represent 12 hour precipitation totals. For the overall correlation, no 
objective limitations are imposed other than the relative humidity threshold. Regional 
constraints are imposed, however, when examining the effectiveness of this technique in 
more localized mountainous regions. 

Perhaps the most obvious limitation of using this quantity to estimate precipitation 
is its exclusion of sub-grid scale (i.e., convective) precipitation processes. While no effort 
is made to quantify the effect of topography on such precipitation, these conditions are 
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Figure 1. Topography (40 km resolution, ın meters) 
C. DATA 
The saturated moisture flux quantity used to represent precipitation is calculated 
from NMC final Global Data Assimilation System (GDAS) analysis fields, which were 


then statistically fitted to hourly precipitation ebsercatienz collected Yo PIOAA, The 
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Figure 3. Synoptically forced 700 mb moisture flux (w0f). 1200 UTC 01 March 
1991. Negative values correspond with upward flux. 
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Figure 4. Total 700 mb moisture flux (wqf), 1200 UTC 01 March 1991. 
Negative values correspond with upward flux. 
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Specific rainfall totals obtained through hourly precipitation data publications 
(NOAA 1991) for the period of study were used to correlate precipitation with the wqf 
values. NOAA’s precipitation totals consist of hourly rain gauge observations made at 
244 locations within Califomia. For our purposes, the observations have been grouped 
in six-hour increments exactly preceding and following the NMC analysis time. This 
gives 12 hour precipitation totals centered on the NMC analysis time which are then 


correlated with tbe corresponding wqf values. 
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IV. SYNOPTIC ANALYSIS 


A. OVERVIEW 

The series of synoptic-scale events that occurred during the period of this study 
were produced by the development of a 500 mb trough off the Northern Califomia Coast. 
The system center migrated northeastward, over the period, and supported a series of 
short-wave troughs and associated surface cyclones (Figures 5 and 6). 

An important factor in the synoptic evolution of this period is the movement of the 
primary jet stream. Approximately midway into the period (1200 UTC 02 March 1991), 
the jet repositioned itself over Central California, from the upper Baja Peninsula where 
it was located early in the period. This shift is evident in the 300 mb isotach analyses 
representative of the jet position at 1200 UTC 28 February 1991 and 0000 UTC 03 March 
1991 (Figure 7). The jet axis shifts from a westerly to a more southwesterly onentation 
coincident with its latitudinal shift. As a result, the weather affecting California was 
dominated by highly baroclinic, fast moving cyclones in zonal flow with the surface 
center well poleward of the main jet early in the period, and weaker, less baroclinic 
systems late in the period, with only a brief transitional time in between. 

As a consequence of the poleward migration of the jet, incoming weather systems 
affected more Northern sections of California late in the period, as evidenced by the 500 


a 


mb height rises observed ever Southein Califereio (tonnes t f eem eom ote most ef the 
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ram Js concentrated jn Southern California early in the period, with virtually no rain 


observed in Southern California by the end of the period. 
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Figure 5. 500 mb heights (meters, left column) and surface pressure, including 
frontal positions (mb, right column) for (from top) 1200 UTC 27 February 1991, 
1200 UTC 28 February 1991 and 1200 UTC 01 March 1991. 
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Figure 6. 500 mb heights (meters, left column) and surface pressure, including 
frontal positions (mb, right column) for (from top) 1200 UTC 02 March 1991, 1200 
UTC 03 March 1991 and 1200 UTC 04 March 1991. 


2 


B. CYCLONE CHARACTERISTICS 


1. 9000 UTC 28 February 1991 - 1200 UTC 02 March 1991 
As mentioned, cyclones affecting the weather in Califomia during this period 
developed on the cold, cyclonic side of the main jet, as seen in the GOES infrared 
imagery at 1200 UTC 28 February 1991 (Figure 8) over the Eastem Pacific and the 
corresponding mean sea level pressure and 1000-500 mb thickness analyses (Figure 9). 


This particular image and analysis are indicative of the general nature of the cyclones 
during this period. The strong surface thickness gradient along 30? N as well as the low- 


level instability upstream of the main frontal cloud bands on the GOES imagery indicate 
the strength of the baroclinic systems passing through California. These systems had 
significant pre-frontal cloud and precipitation bands associated with them. In general, this 
type of pattem ıs similar to that in which Collier (1975) applied his mesoscale 
enhancement model over North Wales, thus providing a good basis for comparison. 
Cyclones of this nature are favorable for observing orographically enhanced 
precipitation, largely due to the forcing being well represented by the large-scale analyses 
and the precipitation primarily associated with the occluded front. However, post-frontal 
convection does pose a problem. While Collier (1975) applied his technique to rain 
associated with both cold and warm frontal passages, Bader and Roach (1977) tried to 
focus only on the mostly stable precipitation produced in the warm sector of cyclones. 
Browning et al. (1974) earlier concluded that orographically induced precipitation was 


most evident m the wann sector, as low-level ee teem plaved ERI ne during and 
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after cold frontal passage. Still, tbew results indicated orographic effects could be 
quantified even during and after the passage of the cold front. Results obtained by Collier 
(1975) showed less degradation of quantifiable orographic enhancement during cold 
frontal passage, as compared to warm sector precipitation, than did the results presented 
by Browning ef al. (1974). Collier (1975) placed more emphasis on distinguishing 
between fast moving baroclinic systems in which mesoscale forcing is minimized, and 
slower, less baroclinic features as far as the observability of orographic precipitation 
enbancement. Since the cyclones in this period of the present study can be characterized 
as occluded surface features with mostly stable precipitation prior to frontal passage, and 
some post-frontal convection. orographic effects should be well represented. This study 
does not examine in detal the effect of mesoscale post-frontal convection on the overall 


correlation between wqf and observed precipitation. 
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Figure 7. 300 mb heights (solid, meters) and isotachs (dashed, knots) for 1200 UTC 
23 February 1991 (top) and 1200 UTC 02 March 1991 (bottom). 
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2. 1200 UTC 02 March 1991 - 1200 UTC 04 March 1991 
The general nature of the three weather systems affecung California during this 
period is represented in Figures 10 and 11. The thickness gradient has shifted northward 
off the California coast and the amplitude of the pre-frontal thickness ridge is less than 


in the earlier cyclones. In addition there is a lack of post-frontal instability depicted on 


the infrared image. An exception to this pattern is the area of cloudiness at approximately 35" 


N, 150? W. As with tbe svstems preceding it, this particular system becomes cut-off 
from the cold air within which 1t developed, and evolves mostly in response to the strong 
southwesterly flow aloft. This system (which reaches California at approximately 1200 
UTC 04 March 1991), like the two seen ahead of it, has little representation of a 
thickness ridge ahead of the front as it reaches land. Without the pre-frontal baroclinic 
forcing, the rainfall associated with this type of system 1s less dominated by stable pre- 
frontal precipitation than cyclones in the earlier period. Due to the uninterrupted supply 
of warm, moist low- and mid-level flow from the southwest during this period, incoming 
weather systems produced less organized and more convective-type precipitation. This 
type of pattern was specifically avoided in terms of quantifying orographically induced 
precipitation by Collier (1975), since the forcing responsible for the distribution of 
precipitation was difficult to represent, even on the mesoscale. Consequently, this period 
of the present study may not be adequately represented by the synoptic-scale and 


orographic effects on precipitation, 
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Figure 8. GOES infrared imagery. 1200 UTC 28 February 199]. 
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3. Summary 


The above synoptic discussion suggests two cases that produced the observed 
rainfall distributions: 
- Case (1) - strongly baroclinic cyclones with a well defined pre-frontal thickness 
ridge. 


- Case (1) - subtropically influenced convective complexes with little or no pre- 
frontal thickness ridge. 


The first case tends to donunate the winter rainy season over California, 
although, the second case 1s not uncommon. Consequently, results from the two periods 


in this study cover contrasting synoptic conditions which produce rain in California. 
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V. RESULTS 


A. OVERALL 


1. Statistics 

Table 1 shows the correlation, slope, and y-intercept from a linear regression 
analysis of wqf and observed 12 hour precipitation for four representative analysis times. 
The correlation values run from a virtually non-correlative 
-0.19 to a marginally respectable -0.55, as shown by the representative scatter plots in 
Figures 12 and 13. The correlation values are grouped at -0.2 and -0.5 for the later and 
earlier synoptic periods, respectively. This suggests that the synoptic and mesoscale 
factors discussed in the previous chapter predetermine the success of this technique to 
be either moderately correlative (-0.5) or essentially non-correlative (-0.2), with little ın 
between. Averaging the slope and y-intercept of all "correlative" times (Table 1 shows 


just representative times) yields -0.010 and +35 respectively. Therefore, 


[-(waf)0.01]+35=P (6) 
equates calculated wqf values with expected 12 hour accumulated precipitation (P), in 
inches. Note that all values of P less than zero represent a precipitation value of zero. 

As will be shown in a later section, correlation is greatly improved in 


mountainous subregions as compared to the state-wide case. Since this technique is 


intended te be “universal” m its application te etve ss elyeetive svneptre -ccale parameters 
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and topograpby, (6) will be appbed in all sub-regions, rather than deriving specific 
statistics for individual sub-regions. 


TABLE 1. LINEAR REGRESSION STATISTICS (WQF VS. 12 HOUR 
PRECIPITATION) 


CORRELATION Y- 
COBEPEICIENT INTERCEPT 
(INCHES) 


! 0000 UTC 28 
FEB 1991 


1200 UTC 01 -.010 +0.35 
MARCH 1991 


0000 UTC 03 #7 -.003 +0.40 
ı MARCH 1991 





1200 UTC 04 -.003 +1.00 
MARCH 1991 


One interesting statistical result is the fairly consistent y-intercept value for the 
linear regression of about +35 (average of six "correlative" times). This indicates that the 
precipitation threshold (the x-intercept) generally corresponds with slightly downward 
saturated moisture flux values. Though the smoothing process tends to distribute upward 
flux values into lee slope regions, it is still likely out-weighed by the "spillover" effect 
observed by Colton (1975). In this, spillover process, upward vertical velocities are 
advected over the summit region, effectively “mixing” upward vertical velocity air into 
lee regions. This process is not accounted for in the wf calculation, which sees all lee 


regions as having dewnward velocities in the absene af avnaptèe colo npward motion. 
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Figure 9. Surface pressure (mb, solid) and 1000-500 mb thickness (meters, 
dashed). 1200 UTC 28 February 1991. 


Additionally, horizontal wind drift, as described by Collier (1975), may be producing 
artificially high precipitation 
observations across summit regions, especially under conditions of strong low-level flow 
with a substantial perpendicular-to-barrier component. 

The other consistent statistical feature is the slope value of either about -0.01 
for the correlated a a or about -0.003 for times showing less correlation. This suggests 


a "flatter" precipitation distribution during non-correlative times, with non-orographic 


processes more dominant SS 
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2. Waf Distribution 
Although the linear regression statistics computed above have a relatively low 
correlation, comparison of horizontal distributions of the wqf-inferred precipitation and 
observed precipitation reveals qualitatively correct distribution. In fact, this comparison 


reveals a significant improvement over NGM model forecast precipitation. 


a. 0000 UTC 28 February 1991 - 1200 UTC 02 March 1991 

Figures 14 and 15 depict the statewide precipitation distributions produced 
by wOf (the synoptic-scale forcing described in Chapter II, NMC's Nested Grid Model 
(NGM) forecast precipitation for 1200 UTC 28 February 1991, actual observations and 
wqf. Of note here is the similarity between the NGM forecast precipitation and the wOf 
field, particularly the maximum over Southern Califomia. This implies that the synoptic 
forcing we are using in our calculations of wgf bears rough equivalency to the NGM 
precipitation. This time is especially informative since the 12 hour 700 mb vertical 
velocity forecast verifies particularly well compared to the kinematically calculated 
analyzed 700 mb vertical velocity (not shown). It should be mentioned that wOf showed 
virtually no correlation with observed precipitation throughout the study. 

In general, we see a significantly greater amount of detail within the 
precipitation field produced by wqf, compared to either wOf or the NGM forecast. Also 
note the "false" maxima produced in the observed field. This 15 mainly a result of the 
multiquadrant contouring technique extrapolating basis functions into regions of sparse 


pman e Mojave Desert, and over the orean) The arepa aside, i tneecture of the wef 
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Figure 10. GOES infrared imagery 1200 UTC 02 March 1991. 


distribution reflects the observations reasonably well, particularly the maxima observed 
over the mountains of Southern California and the Sierra Nevada. Note that correlations 
are performed between the gridded wqf fields and station locations (rather than gridded 
precipitation). It should be noted that overall correlation is generally constant from 0000 
UTC 28 February 1991 to 0000 UTC 01 March 1991, with the exception of time 0000 
UTC 01 March 1991, which has a correlation coefficient of -0.21. This drop in 


correlation 1s most likely in response to the low-level convection seen in Figure 8 located 


between the main frontal bands. This pattern ar omo iue pow ane IS n 


12 


post-frontal precipitation distribution ıs determined less by orography than convective 


processes. 


b. 1200 UTC 02 March 1991 - 1200 04 March 1991 

Figures 16 and 17 depict the same parameters as Figures 14 and 15 
except for 1200 UTC 04 March 1991, which is representative of the second period. Here 
wqf is generally a less accurate representation of observed precipitation, particularly along 
the Central California Coast. Stil, orographic maxima produced by wgf over the 
mountains of Northem Califomia and the Northern Sierra Nevada are reasonably 
representative of observed maxima. Though the corresponding NGM precipitation 
forecast shows good areal coverage, 1) 15 devoid of any of the above structure. In this 
case the NGM vertical velocity forecast. on which precipitation 1s largely based, shows 
a strong area of upward motion Jocated over Northern California, whereas the verifying 
kinematically calculated analysis shows this area of upward motion to remain off-shore 
(not shown). As such, the NGM precipitation field owes its accuracy more to chance than 
to improved dynamical representation. It seems that neither the NGM forecast or the 
GDAS analysis fields axe able to pick up on the forcing responsible for the locally heavy 
rainfall in this period. Referring back to Figure 6, we see an analyzed ‘front’ off the 
Northern California coast at this time, with most of the state within the warm sector. Á 
random check of station observations supports the strong southerly flow depicted by the 
GDAS analyses. Convective lifting would. in this case. be the best explanation for the 


Sup ade AN nl liste los ne 8: As nete: 
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Figure 11. Surface pressure (mb, solid) and 1000-500 mb thickness (meters, 
dashed). 1200 UTC 02 March 1991. 


earlier, this period is dominated by only marginally baroclinic features in the moist 
southwesterly flow aloft. As such. convection is likely to be equally important in both 
the pre- and post-frontal periods. Though convective processes are parameterized in the 


NGM precipitation forecasts. they are not accounted for at all in the calculation of wef. 
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3. Correlation Limiting Factors 


a.  Synoptic Conditions 

As described by Collier (1975), synoptic conditions favorable for 
observing orographic enhancement of precipitation are characterized by fast moving 
baroclinic systems. The thinking is that slower moving, or non-baroclinic systems allow 
more non-orographic and/or mesoscale mechanisms to influence rainfall distribution. This 
is not to say that topography does not influence mesoscale structure. Rather, since only 
the large-scale structure is considered. topography is only allowed to influence the 
synoptic-scale fields. As mentioned previously, the period prior to 1200 UTC 02 March 
1991 was characterized by baroclinic features well represented by the synoptic fields, 
whereas non-orographic and/or a mesoscale processes dominated after this tune. The 
Statistical correlations support this characterization of the dominant precipitation type 
during the perieds. However, as seen in the horizontal distribution of wqf during the 
second period, the qualitative distribution of precipitation is realistic in the second period. 
This suggests that the amount of precipitation is poorly represented in more convective- 


type periods, but the general distribution is still strongly influenced by the topography. 


b. Smoothing 


Aside from subgrid-scale effects, smoothing of the topographic w used 


to calculate (wg), seems to be most detrimental to correlation values. As previously 


uoussedeimiesoseste orecssphic feremesrocduee In esse sinn Scnepric-scale 1000 


mb winds over 4Q km resolution tepographiy refes tio-pags semoethine to be 
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Figure 12. Observed precipitation (vertical axis, inches) vs. wqf (horizontal axis) 
and calculated linear regression for 0000 UTC 28 February 1991 (top) and 1200 
UTC 01 March 1991 (bottom). 
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compatible with the synoptic-scale forcing. As a result, stations upwind of topographic 
barriers, such as Fresno (station FRES, Figure 18) are generally assigned wqf values far 
in excess of the actual calculated forcing. Note (Figure 18) that the smoothed waf value 


for FRES (-15) is nearly twice that of the unsmoothed value (-8). 
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c. Location 

Certain Jocations, such as coastal areas, are more likely to have higher 
actual relative humidities than inland areas, which are analyzed as equal on the synoptic- 
scale. As a result, more of the vertical moisture flux waf is likely to be converted to 
precipitation in coastal areas than inland areas, where a larger amount of moisture flux 
(wqf) may be required to maintain or reach saturation. This relationship is evident in 
comparing the linear regression corresponding to the San Bernardino - San Gabriel 
subregion (subregions are defined in the following sections) to that of the Southern Sierra 
Nevada (Figure 19). Note that a wqgf value of -100 represents approximately a 2.75" 12 
hour precipitation accumulation according to the linear regression derived for the San 
Bemardino/San Gabriel subregion, but less than 1.0" based on the Southern Sierra Nevada 


subregion’s statistics. 


d. Other Considerations 
Obviously, numerous other considerations affect the correlation of 
precipitation and 700 mb saturated moisture flux. Redistribution of moisture above 700 
mb, quantification of evaporation, and wind channeling induced by the topography are a 
few examples. Though important, such processes cannot be easily accounted for on the 
synoptic scale. In fact, these factors will require more specific consideration if wqf 
cannot adequately represent the forcing patterns which lead to precipitation distributions 


in mountainous regions. 
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Figure 13. Observed precipitation (vertical axis, inches) vs. wqf (horizontal axis) 
and calculated linear regression for 0000 UTC 03 March 1991 (top) and 1200 UTC 


04 March 199] (bottom). 
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Figure 14. NGM 12 hour precipitation forecast, valid 1200 UTC 28 February 
199] (mm, top) and wOf (bottom) fox 1200 UTC 28 February 1991. 
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Figure 15. Wqf for 1200 UTC 28 February 1991 and observed precipitation 
(inches, bottom), 0600 UTC 28 February 1991 - 1800 28 February 199]. 
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Figure 16. NGM 12 hour precipitation forecast valid 1200 UTC 04 March 91 
(mm, top) and w0f (bottom) for 1200 UTC 04 March 1991. 
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Figure 17. Wof (top) for 1200 UTC 04 March 199] and observed precipitation 
(inches, bottom), 0600 UTC 04 March 1991 - 1800 UTC 04 March 1991. 
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Figure 18. Unsmoothed (top) and smoothed (bottom) orographically forced 
vertical velocity (cm/s). 0000 UTC 01 March 1991. 


B. REGIONAL 


J. Approach 

While the general results are indicative of the overall value of the wqf 
parameter in representing precipitation over a wide region, it is perhaps more important 
to determine the extent to which the correlation can be improved when applied to a 
specific topographic region. In this way, wqf may be able to identify specific regions or 
even drainage basins likely to receive heavy precipitation under certain forecast 
conditions. For this reason, the specific perfonnance of wgf in regions of well defined 
topography is examined to determine the ability of wgf to represent precipitation in 
specific regions. In effect, the extent to which synoptically derived low-level winds 
interacting with topography can represent the mesoscale distribution of precipitation 
within mountainous regions will be determined. Little or no improvement of the regional 
correlation coefficient over the general case would suggest limited usefulness of the 
synoptic-scale winds in determining orographic forcing. An increase in correlation would, 
however, indicate that a significant component of the precipitation forcing is due to the 
large-scale wind field interacting with regional topography, and that this quantity is more 
easily isolated in regions of well defined topography. Specific correlation coefficients are 
presented as required. in each subregion. 

The regions chosen for this specific purpose all have topography which is 
relatively well represented at 40 km resolution. In each case. regions are defined by the 


E CUI (ODOT TN keatines heise salat nat es I >. tete Results are 
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presented by showing tbe observed precipitation across a topographic cross-section which 
is Chosen based on observational density as well as topographic gradient. Observations 
are then compared to the wqf derived quantities. Topography requiring finer resolution to 
be represented does not fare well, as smoothing tends to destroy the forcing maxima (a 
good example is the Santa Cruz mountains). As much as possible, limitations, as well 


as strengths, corresponding to specific regional correlation will be identified. 
2. Topographic Subregions 


a. San BernardinolSan Gabriel 

This subregion was particularly hard hit by storms early in the period. 
Figure 20 shows the regional topography to consist of a broad flat area bordered to the 
northeast by the abrupt rise of tbe San Bernardino range. Figure 21 shows the actual 
observed precipitation accumulations along with the wqf derived quantities, for 0600 UTC 
-1800 UTC 01 March 1991 combined with a representative cross-section of the regional 
topography. This particular cross-section is representative of the general correlation 
within this subregion, and was chosen since it represents the period of maximum rainfall. 

The wqf parameter seemed to do a particularly good job handling the 
Mojave Desert stations throughout the period, while showing only moderate success in 
the upslope regions. This may best be explained by the coarse resolution of the 


topography southwest of the crest. Our technique sees only a broad, continuously sloping 
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Figure 19.. Waf vs. observed precipitation (inches, vertical axis) for San 
Bemardino/San Gabriel subregion (0000 UTC 28 February 1991, top) and Southem 
Sierra Nevada subregion (1200 UTC 01 March 1991, bottom). 
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region, whereas many reporting stations lie mn dissected valleys (i.e. Riverside at 301 m 
is assigned an elevation of 550 m by the 40 km topography). Also note (stations 1, 2 and 
3 on Figure 21) the inflated wgf values close to the ocean. This results from the 
smoothing process described earlier. 

With some notable exceptions (1.e., Lytle Creek Ranger Station (LYCR 
#11)), wgf does a decent job of representing the regional precipitation showing a 10% 
increase in the correlation coefficient over the state-wide case. Forcing responsible for 
the anomalous observations (i.e., LYCR and CRES) is certainly of too small a scale to 
be accurately accounted for by this technique. Stil, the ability of wqf to correctly 


1 


produce a specific area of heavy precipitation (21.5" in 12 hours), as well as define the 
lesser amounts observed in the Mojave Desert (stations VICT, BORO, DAGG and MOJA) 
is encouraging. By comparison, NGM 12 hour precipitation totals forecast for this region 
wete all less than 0.4 inches with no upslope or lee structure in the NGM precipitation 
(not shown). 

Of course, wqf has the benefit of being statistically fitted to actual 
observations, where the NGM forecast represents a strictly physical quantity. 
Nonetheless, the statistical fit achieved by waf is derived objectively, and is intended to 


be applied in all conditions. Still, an application of this method to independent data is 


needed to determine predictive skull. 
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Figure 20. Topography (meters) and station names for San Bernardino/San 
Gabriel subregion. Precipitation amounts (inches) are for 12 hours ending 1800 UTC 
01 March 1991. 
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Figure 21. Topography (solid line, meters), wqf precipitation (dashed, inches), 
and observed precipitation (points. inches) for San Bemardino/San Gabriel subregion 
(12 hour period ending 1800 UTC 01 March 1991). 


It should be mentioned that the significant rain events in both this 


subregion and the Santa Barbara subregion which follows ended by 1800 UTC 02 March 


1991. As a result, the conditions which lead to deterioration of correlation values late in 


the period were not observed in these regions. 


b. 


good waf/precipitation conslation «Figmae., 


Santa Barbara 


The Santa Barbara subregion (Figure 22) 1s characterized by particularly 


menn 


50 


N recat ys ' eve cus Yi tan nose fr. 


..f 


: ranging from 


-0.52 to -0.85 for the different analysis times. Note that correlation during the post- 
frontal convective period of 1200 UTC 28 February 1991 - 0000 UTC 01 March 1991 
shows less deterioration in this region than ın the state-wide case (-0.52 vs. -0.22). 
Simuar to the previous subregion, lee stations (NE) are represented highly accurately, 
with the majority of the variance seen in the upslope regions (SW). Again, subgrid-scale 


topography likely produces large vanations over small distances. 
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Figure 22. Topography (meters) and station names for Santa Barbara subregion. 
Precipitation amounts (inches) are for 12 hours ending 1800 UTC 01 March 1991. 


A good example of an error caused by insufficient resolution of the 
topography is evident at time 0000 UTC 28 February 1991. Figure 24 shows the wqf 
maximum to be "transposed" to the southeast from the observed maximum, therefore 


limiting any correlation in thia location, Interestinzh the combed ost topegraphy in tus 
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Figure 23. Topography (solid line, meters), wqf precipitation (dashed, inches), 
and observed precipitation (points, inches) for Santa Barbara subregion (12 hour 
period ending 1800 UTC 01 March 1991. 


region (not shown) depicts a narrow northeast trending ridge in the vicinity of Thousand 
Oaks (station LECH), northwest of which is the Ventura Valley (station SIMI) This 
topography is represented as a continuous slope at 40 km resolution, however. As such, 
the strongest orographically forced component of waf is located over the Ventura Valley 
(near SIMI), whereas in reality this forcing maximum is located between the Ventura 


Valley and Pine Mountain (station PINM). As the southeasterly winds become more 


southerly after this fime. wtf and olo rrr oT Seggans Y 5 wee rry AZI e + آ5‎ much POTE 


53 


a 
- 


1 
' 
, 
1 
1 
۷ 
0 





Figure 24. Wqf and observing stations with 12 hour precipitation totals ending 
0600 UTC 28 February 1991 (inches) for the Santa Barbara subregion . 


favorably, as the large-scale winds interact with topography that is better represented at 
the 40 km resolution. 


Overestimation of precipitation at the coast (Figure 23, stations 1-5) does 
not seem to be as premalent mts cio" 


Bar o d LS EU eye el, 


of the extreme 
slope at the Jand-ocean margin, Specifically, although omeothing of the topographic w- 
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field stil! produces upward orographic forcing well into oceanic regions, the rapid increase 
in elevation onshore produces orographically enhanced observed precipitation amounts 
even near the ocean. Moreover, the apparent high coastal correlation is more likely a 
coincidence as opposed to an increase in skill over, for example, the San Bemardino 
region. As with the previous sub-region, no significant precipitation was observed after 


1800 UTC 02 March 1991. 


c. Southern Sierra Nevada 
The most notable feature seen in this region (Figure 25), commonly 
koown as the High Sierra. is the Jack of observations. Still we can make at least one 


important observation: 
- Wof generally overestimates precipitation for these inland locations. 


This trend is evident throughout the study period and is represented in 
Figure 26. We can generally note that wqf responsible for a specific ram amount along 
the coast (ie., Santa Barbara) represents a somewhat lesser amount at more inland 
locations. This feature is almost certainly a function of the available moisture, such that 
inland locations require a greater wqf to maintain saturation than is required near the 
coast. No effort is made to quantify this relationship, however it must be considered 
when using wqf based on inter-regional statistics (similar to adjusting Model Output 


Statistic (MOS) precipitation forecasts for a specific location). 
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Figure 25. Topography (meters) and station names for Southem Sierra Nevada 
subregion with 12 hour precipitation totals (inches) ending 1800 UTC 01 March 
199 


Note that Figure 26 shows wqgf and observed precipitation both before and 


after | 200 PDC: OZ Macb A Se pe ua, ا‎ res a n nan FO on me st RZ See 


0.70 respectively, Thongh in the latter, a large regres of correlation ne mans, nep- 
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Figure 26. Same as in Figure 23 except for Southern Sierra Nevada subregion 
(12 hour period ending 1800 UTC 01 March 1991, top, and 0600 UTC 04 March 
1991, bottom). 
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otographic effects are more evident, especially considering the rough equivalency of 
precipitation totals at the San Joaqum Valley station Corcoran (CORC) and amounts 
recorded higher m the Sierra. “The low slope value (-0.004) of the linear regression 


supports this assignment. 


d. Northern Sierra Nevada 

This subregion (Figure 27) depicts tbe gradually sloping region of the 
Sierra Nevada between Sacramento and the Lake Tahoe basin. The chosen cross-section 
(Figure 28) coincides roughly with the American River drainage, which is represented 
particularly well by the 40 km topographic resolution. As a result, we see very good 
waf/precipitation correlation in this sub-region prior to 1200 UTC 02 March 1991, with 
correlation coefficients between -0.70 and -0.80. As before. we note the overestimation 
ef precipitation in the Sacramento valley (WHEA), as well as a general tendency to 
underestimate precipitation (slope values average -0.006 in this subregion). The generally 
good fit seen in the first period at the summit and lee stations of Soda Springs (SODA) 
and Truckee (TRUC) respectively represent the ability of this technique to account for lee 
“spillover”, due to the statistical correlation producing a positive y-intercept. 

Perhaps the best example of subgrid-scale processes affecting the 
distribution of precipitation is seen at time 1200 UTC 04 March 1991 (Figure 27). Note 
the extracrdinarily high precipitation observations over Donner Summit, including 
Truckee, as compared to precipitation estimated by wgf (regional correlation coefficient 
adi slope values aie =O ane e ML LL fe we kan thee so sp um cu 


resolving the wind and moisture field anomalisg resgpyonsible for tbe: deztribytiopn, This 
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example brings into question wgqf's level of usefulness given this type of synoptic 
situation. 

Of particular note concerning this region is the ability to compare results 
obtained here to those obtained over a similar region by Colton’s (1976) 4.3 km model. 
Though Colton’s (1976) use of an hourly precipitation rate prohibits direct comparison 
to the 12 hour rate presented here, we can make some qualitative assessments. Other than 
Colton (1976) showing considerably greater detail within the upslope region, including 
the secendary maximum produced by Mt. Rose rising to the east of Lake Tahoe, wqf 
produces a similar enhancement pattern, prior to 1200 UTC 02 March 1991, in terms of 
tbe location of maxunum as well as the intensity. Of further note is that even at the 4.3 
km grid scale, Colton’s (1976) calculations at lee stations generally underestimate 
precipitation, even with a "spillover" term included. However, Colton’s (1976) fine mesh 
resolution allows a better representation of convective forcing. In these types of 


situations, wdqf is far inferior to Colton’s (1976) model. 
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Figure 27. Topography (meters) and station names for Northern Sierra Nevada 


subregion with 12 hour precipitation totals (inches) ending 1800 UTC 04 March 
l 


e. — Shasta/Siskivon 
The final sub-region (Figure 29) contains not only the most complicated 
topography. but also wag subiect te the svideost cose Co lC STEM 


Figure 39 shows representative waf/precipitation values before and atte: 1200 UTC 
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Figure 28. Same as in Figure 23 except for Northem Sierra Nevada subregion 
(12 hour accumulations ending 0600 UTC 02 March 1991 , top, and 1800 UTC 04 
March 1991, bottom). 
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02 March 1991 with corresponding correlation coefficients of -0.68 and -0.46 respectively. 
Also, slope values for this subregion average -0.012, which, being relatively high, is 
consistent with its more coastal location. Despite the obvious smoothing of topographic 
features at this resolution, the overall wqf representation of precipitation is quite 
encouraging, in particular the secondary maximum seen between Happy Camp (HAPP) 
and Elk Valley (ELK V). Even more interesting is the less obvious loss of correlation late 
w the period. This may be due to this subregion’s northemmost location. Here, 
synoptically derived subsidence 1s less than that at more southern locations at this time, 
allowing wqf to "create" upward motion by offsetting the synoptic-scale subsidence. 
Clearly, the improved correlation is only coincidental since convective processes are not 
accounted for in any region. Another possibility is more consistent forcing. The 
proximity of this subregion to the surface low pressure systems late in the period may 
have resulted in a more consistent low-level southerly flow. In this situation, though 
convection is still unaccounted for, the synoptic- scale winds are probably acceptable 
representations of the low-level flow in this subregion. This, combined with the strong 
tendency for the north end of the Sacramento Valley to channel winds with a southerly 
component may have allowed wgf to reach a level at least of the same order of magnitude 
as the convective forcing. As such, this period’s correlation is likely more an effect of 
the particular local situation interacting with the synoptic scale, rather than wqf being 


more accurate in this subregion 
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Figure 29. Topography (meters) and station e 
ipitation totals 


names for Shasta-Siskyou subregion, 
with 12 hour precip (inches) ending 1800 UTC 01 March 1991. 
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Figure 30. Same as in Figure 23 except for Shasta-Siskyou subregion (12 hour 
accumulations ending 1300 UTC 01 March 1991, top. and 1300 UTC 04 March 
1991. bottom). 
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f.  Non-mountoinous Regions 

Though this study is aimed at representing mesoscale precipitation within 
mountainous regions, it is also hoped that including the effects of local topography will 
improve distribution in low lying regions, or regions in which topography is sub-grid 
scale such as the Santa Cruz Mountains (Figure 31). Without going into great detail, we 
can qualitatively assess the success of wqf as weaker where topography is not as well 
defined. As before, regions of less dramatic topography are most susceptible to sub- 
synoptic processes dominating the local distribution of precipitation. This is especially 
evident when surface forcing is not represented well by the synoptic analysis (i.e. upper 
level features and convection). 

As far as mountamous regions which are poorly represented at 40 km 
resolution, there is sufficient evidence to suggest improving topographic resolution in 
order to identify waf forcing in these regions. Again we refer to the precipitation 
distribution within the Santa Cruz mountains (Figure 32). Present resolution does not see 
the steep, narrow nature of this range, but observations support strong local orographic 
influence. Increasing topographic resolution may provide sufficiently large quantities of 
u'gf to "survive" the smoothing process realizing that this representation will only be of 
use if the synoptic analysis is at least close to reality, which is certainly far from 


guaranteed. 
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Figure 31. Topography (meters) and station names in the vicinity of the Santa 
Cruz Mountains with 12 hour precipitation totals (inches) ending 1800 UTC 28 
February 1991. 


3. Comparison to NGM Precipitation Forecasts 
Easlier we established the rough equivalency of the NGM precipitation forecast 
with the term used bere to represent synoptic-scale forcing, wOf. We have seen examples 
of the improvement the topographic forcing term has on the wf, field in representing 
observed precipitation during the first period. “This is essentially the same level of 
improvement wgf wall hielo shese eyer the IAS ld AOS allowing 


for the NGM convective parameter (thig parameter could juzé 23 eat lv be bout 
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Figure 32. Wgf and observing stations with 12 hour precipitation totals ending 
1800 28 February 1991 (inches) in vicinity of Santa Cruz Mountains. 


into wgf). Note that while the NGM and w0f are able to define general precipitation areas 
reasonably well, there is virtually no structure other than that associated with tbe synoptic- 
scale forcing maxima (usually coincident with 700 mb vertical velocity), and therefore 
weak correlation with observations. 

Late in tbe period. NGM represents precipitation with only slightly less overall 
accuracy than wqf. Recall that for these times, the NGM vertical velocity forecast 
predicted upward motion in areas which verified to be subsident. Still. its convective 


parameter may Bere ent Reel rn trees 
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VE. CONCLUSIONS 


A. APPLICATION TO FORECASTING 

Tt is well known that precipitation amounts can vary dramatically over a small area. 
Inclusion of topography into fine-mesh models has shown favorable representation of this 
variability (Colton 1976, Collier 1975). Unfortunately, the data network needed to run 
tbese models is not generally available to the forecaster, who can only combine local 
knowledge with the svnoptic-scale products he receives. In this thesis we have attempted 
to provide a quantitative estimate of mesoscale precipitation distributions over 
mountainous terrain, without requiring mesoscale information, other than topography. 
Conclusions can be summarized by the following points: 

- Given that synoptically induced forcing is well represented in the analyzed or 
forecast fields, wqf can predict a substantial degree of the mesoscale precipitation 
distribution. 

- Most of the discrepancies between wqf and the observations occur within upslope 
regions (as opposed to lee stations) due to unresolved small-scale topographic 
forcing. 


- Correlation is greatly unproved when considering well defined topographic regions. 


* Wdqf shows significant unprovement over the NGM precipitation forecast fields. 


Coincident with our conclusion that wqf is a good indicator of precipitation under 


well defined svnoptic-3cal- kesye tre Gesprelettesro: eo TE ee een n Jane ox rs anf less value 
when non-oregraphic and/or mesoscale piersia deminate precipitation forcing. 
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Therefore, to be of use to the forecaster, he/she must first determine the applicability of 
waf based on the synoptic situation. There is sufficient evidence to suggest wqf may be 
of some value even when no synoptic-scale forcing is indicated. In this case, the 
forecaster must weight the waf field accordingly, but be aware that areas of consistently 
strong flow may still be represented by wqf precipitation maxima. As such, assuming wqf 
is produced as a forecast field, the forecaster must identify its credibility based on the 


following parameters: 


- Presence of a well-defined thickness ridge ahead of a front. 


- The potential for unstable precipitation. 


Jf these conditions can be correctly evaluated, wqf will be sufficiently useful to the 
forecaster ın predicting locally heavy precipitation to warrant its calculation by regional 


forecast offices. 


B. RECOMMENDATIONS 

The most immediate recommendation would be to incorporate actual NGM forecast 
fields into the wgf calculation. This would provide a basis for verifying the wqf 
precipitation against the actual NGM forecast precipitation amounts. Such an application 
will be able to more definitively evaluate wqf as a forecast parameter. 

Secondly, there may be some benefit to computing wqf over even finer topographic 
resolution (1.e.. 20 01 even 10 kpo gines the foster e bi b ect relation will be 


prevalent regardless of the topographic resclutien Such an mersa + ° resolulion mav 


improve wqf representation of precipitation by pinpointing areas of maximum 
enhancement and by reflecting the orographic influence of steep topography with less 
borizontal extent. There will, however, most certainly be a point at which no increase in 
wqf accuracy is seen with increasing topographic resolution. Furthermore, objective 
restrictions potentially can be applied to the calculation of wqf in flat regions which are 
adjacent to significant topographic barners. The intent here is to limit the forcing to the 
actual topographic slope rather than distributing it over a larger region as done by the 
current smoothing process. 

As a final note, I wish to empbasize that no new ground has been broken here as 
far as our understanding of the precipitation process. I bave simply adapted seeder-feeder 
theory of orographically enhanced precipitation to a format which only requires synoptic 
scale information. Though the results are far from ideal, I think the technique exhibited 


sufficient success to warrant jts consideration by the forecaster. 
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